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Array detector-based instruments are now fundamental to measurements of ozone and 
other atmospheric trace gases from space in the ultraviolet, visible, and infrared. The 
present generation of such instruments suffers, to a greater or lesser degree, from 
undersampling of the spectra, leading to difficulties in the analysis of atmospheric 
radiances. We provide extended analysis of the undersampling suffered by modem satellite 
spectrometers, which include GOME, SCIAMACHY, OMI, and OMPS. The analysis 
includes basic undersampling, the effects of binning into separate detector pixels, and the 
application of high-resolution Fraunhofer spectral data to correct for undersampling in 
many useful cases. 
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1. Introduction 

Array-based spectrometers used in atmospheric remote sensing can suffer substantially from 

spectral undersampling, with negative consequences to the quality of data retrieved from the 
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measurements. This was first obvious in measurements from the Global Ozone Monitoring 
Experiment (GOME), 1 where scientists fitting trace gases from GOME spectra found quite large 
systematic fitting residuals, and large fitting errors, even after correcting for Doppler shifts 
between radiances and irradiances. Ref. 2 recognized that this was mainly due to spectral 
undersampling and presented a technique for correcting most (> 90 %) of the undersampling error 
in spectra] regions where atmospheric absorption effects are small. For GOME, this includes 
fitting regions used for nitrogen dioxide (NO2), bromine monoxide (BrO), chlorine dioxide 
(OCIO), and formaldehyde (HCHO). If used with caution, the technique can also be applied 
successfully to O3 and SO2. (Caution is required because the absorption optical depths for O3 are 
higher than for the other trace species listed; correction assumes that, to first order, the radiance 
spectrum consists mostly of back scattered Fraunhofer structure. The SO2 absorption occurs in 
the region where O3 also absorbs strongly.) The technique consists of comparing fully-sampled 
and undersampled versions of a high resolution Fraunhofer reference spectrum, 3 with the 
difference being the effect of undersampling. Ref. 4 suggested that the observed residuals are 
induced by the resampling required to compare Earth radiance and solar irradiance spectra in the 
fitting, because they are measured with different Doppler shifts of the ERS -2 satellite with 
respect to the sun. It is demonstrated here that this is not entirely the case: Wavelength shifts 
between GOME radiances and irradiances are larger than can be accounted for by Doppler shifts. 
Ref 4 also presented a version of the technique of Ref. 2 implemented in the Scanning Imaging 
Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) operational processor. 
The technique is now also used for measurement of HCHO 5 and NO2, 6 and has been 
implemented in the operational processor for measurements of BrO, OCIO, and HCHO by the 
Ozone Monitoring Instrument (OMI). 7 
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The present analysis improves on the previous understanding of correction for 
undersampling by fully considering the sampling theorem in conjunction with the instrument slit 
function. The analysis is applied to measurements by GOME and OMI, since they represent the 
two major instrument types: diode array detectors (GOME) and CCD detectors (OMI). The 
development applies as well to other instruments, including SCIAMACHY (diode array 
detectors) and the Ozone Mapping and Profiler Suite (OMPS), part of the National Polar 
Orbiting Environmental Satellite System (NPOESS), which uses CCD detectors. 

2. Definitions and Assumptions 

The spectrometers considered here use array detectors, which respond to the incoming light over 
finite regions, with defined response profiles. We will use the term instrument line shape (ILS) 
for the response of a spectrometer to a monochromatic source up to where it enters the array 
detector, and the term instrument transfer function (ITF) for the ILS convolved with the detector 
pixel response. The following properties are assumed to be valid over the region of the array 
detector that is necessary to consider for ITFs and undersampling corrections at a given pixel, as 
developed here. 

1. Linear response. Each pixel responds linearly to input light intensity. 

2. Equal pixel response. Pixels respond equally to photons of different wavelength. Pixels 
respond (or can be calibrated to respond) equally in output signal for equal input light 
intensity. 

3. Linear dispersion. Pixel wavelengths are at equal increments. 

4. No endpoint issues. The region under consideration is sufficiently far from the array end 
for wavelengths beyond the array wavelength range to contribute to the signal. 
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(Proximity to the array end would require re-calculation of the ILS from deconvolution of 
the ITF, as discussed below.) 

5. Calibration light source issues. Instrument transfer functions (ITFs) are determined 
using either 

• Lines with negligible spectral width or; 

• A tunable source with negligible spectral width and equal intensity as it is tuned 
over the ITF. 

6. Slit width variation is negligible over the ITF centered at each particular wavelength. 
Higher-order corrections will likely be needed in the future to account for the breakdown of 
some of these assumptions. 


3. Spectral Undersampling 

Consider a wavelength range extending from X min to )i max that is fully or partly sampled by an 
array of detectors, spaced at wavelength increment 6/ 1. Any continuous incoming signal spectrum 
extending over the wavelength range can be represented fully by expansion in spatial (i.e., 
wavelength) frequency in a Fourier series: 


S(A) = a 0 + Eh. cos kt + b k sin kt] ~a 0 + X K cos (a) k (A - A mm )) + b k sin(<y t ( A - A mm ))], 

k~\ k= 1 
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and the spatial frequencies are: 
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Suppose that the signal is band-limited , that is, limited in content of spatial frequency 
components to a maximum spatial frequency co, mu . Then, by the sampling theorem, ' the 
information content of the spectrum is fully known if the spectrum is sampled over the full range 
A„„„ to X max to twice this maximum spatial frequency, 2 a> max (the Nyquist sampling frequency)-. 


SA 
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In this case, the signal expansion includes only the terms necessary to measure spatial 


frequencies < (o inax 


N 

S(A) = a 0 + 2k coskt + b k sin fa], N = 

k = I 
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The spectrum is completely determined by this expansion. Its values at points other than the 
sampled points can be determined using the fact that the value at each sampled point X n 

Q 

represents the intensity c„ of a sampling function, sine [2k (X n - A) / <5A], centered at that point. 


sme 


[2n(A n -A)/SA] = sin fW n 
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plus the constant offset ao : 

C. = SO,y - a, = X K cos kl. + b t sin k,, ], t, = ZZ&zi & l 
k=\ A/l 


(7) 


The sampling function is written here in this way to emphasize the pixel dependence and 
orientation of the following discussion. The sampling function presented here is actually an 
approximation to a fuller and more complex form, 9 ’ 10 but the difference amounts to a completely 
negligible correction except within several sample points of the A m ,„ and X max endpoints. 

If the signal is not band-limited to co max but is still sampled to only 2co max , spatial 
frequencies greater than com^ are aliased into the band 0 < co < co max , with signal information for 
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gw < (o < 2 (Omax appearing at ow - co, information for 2&W < to < 3to max appearing at co - 
2 co, no*, etc. (Ref. 9, pp. 16-18). It will be demonstrated later that information aliased from the 
spatial frequency band ow < co < 2cOmax represents the most important sources of problems for 
GOME, and the main target of corrections. Problems from aliasing in the spectra become most 
evident when it becomes necessary to resample a spectrum in wavelength, for example when 
comparing an atmospheric radiance spectrum with a solar irradiance spectrum to determine 
atmospheric composition from molecular absorption lines, but the interference from aliasing is 
present in any case when the spectrum is undersampled (i.e., the spectrum is not fully Nyquist 
sampled). For example, synthetic spectra calculated during the fitting process to determine 
abundances of atmospheric gases would not normally include aliasing, while the measured 
spectra would. 

A Nyquist-sampled spectrum 5 (A) is fully described by summing over the contributions 
from the m individual sample points: 

SU) = c 0 + J c, sinc[2*(4 - A) / SA). (8) 

i=i 

It is now possible to investigate the case where the spectrum is not fully Nyquist sampled. 
Consider a spectrum input to the instrument S mp (A.). If the spectrum is completely known a 
priori , it can be expanded as before in a Fourier series 

cos kt + b k sin kt\ (9) 

*= i 

Si„p can be separated into a Nyquist-sampled part, S„ yq , containing only spatial frequencies < 
comax, and an undersampled part, S un d'- 

N 

S w (/l) = ao+Xk cos kt + b k sin kt\ (10a) 

t=i 
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(10b) 


S tmd U) = £ [a t co&kt + b k sinkt\ 

k = N + \ 

The undersampled part of the spectrum is thus 

m 

S un M) = V a >-£ c < sinc[2^a, - A)ISA.)-c 0 , 

1=1 

( 11 ) 

where m runs over the sampled points. 

4. Slit Functions (Instrument Transfer Functions) 

An ITF serves as a low-pass filter to limit the spatial frequency content of the spectrum. Ideally, 

an ITF would limit the spectral information to frequencies < co max . The ITFs for satellite-based 

spectrometers are measured in one of three manners: 

1. By using a reference line lamp (often a PtNeCr lamp 11 ). This method has the advantage 
that lamp lines are much narrower than the ITF. The disadvantages are that lines are not 
always completely separated, that spectral coverage may be inadequate in some regions, 
and that a set of points (one per detector pixel) is mapped out, rather than a continuous 
ITF. This method was used for GOME and SCIAMACHY. 

2. By using a tunable source consisting of a broadband light source and a monochromator. 
This method has the advantage of mapping out a continuous ITF, but the disadvantage of 
being a spectrally broader source than the PtNeCr line source. This method was used for 
OMI, where the full-width at half-maximum (FWHM) of the source was < 0.1 OMI 
detector pixel (M. Dobber and R. Dirksen, private communication, 2003). The residual 
effects from finite source width could, in principle, be reduced using the deconvolution 
techniques introduced below. 
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3. By fitting flight spectral irradiance data to a high-resolution solar reference spectrum, 3 
where the fitting includes wavelength adjustment and simultaneous fitting to a 
parameterized 1TK 2 ' 12 

An ITF, T (A), can be expanded to include Nyquist-sampled and undersampled portions, ignoring 
(for now) detector pixel binning: 

HA) = y 0 + J>, sinc[2^(/l ( -X)/SX]+ T^X). (12) 

i— #i 

The limits a and b are selected to include portions of the detector array where the slit function 
contributes significantly. As examples, we show the Gaussian GOME ITF we normally use for 
wavelength calibration purposes. A more complex, compound hypergeometric, ITF was 
determined during the instrument characterization. 13 We find that the Gaussian spectrum 
provides better wavelength calibration for GOME spectra, and is used routinely in our data 
analyses. It also provides much better undersampling correction for GOME spectral fitting. We 
show an OMI ITF (M. Dobber and R. Dirksen, private communication, 2003) for this same 
wavelength region, and for the NO 2 fitting region (405-465 nm; the ITF determined at 432 nm is 
used here). The OMI ITFs are selected for CCD row 150, corresponding to a viewing azimuth 
angle of 29.5°, about half way from the nadir view to the extremity of the OMI swath. Additional 
types of ITFs for OMI may be considered once there is flight data for comparison. For all ITFs 
shown in this section, the complication of binning over the detector pixel response function is not 
yet included. This will be discussed in a later section. 

Figure 1 shows the GOME Gaussian ITF, determined from ERS-2 orbit 81003031 
(October 3, 1998), used here as a test orbit for determination of undersampling correction, 
appropriate to the fitting window used for BiO retrieval in GOME. 2 The FWHM of the ITF is 
0.160 nm, and there are 1.4 samples per FWHM in this region of GOME spectra. The ITF is 
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additionally decomposed into fully-sampled and undersampled components. Figure 2 shows this 
ITF with the hypothetical sampling to twice the GOME sampling frequency. Figure 3 shows the 
OMI ITF for this same wavelength region (FWHM = 0.421 nm, 2.8 samples per FWHM), and its 
decomposition into fully-sampled and undersampled components. Figure 4 shows this ITF with 
the hypothetical sampling to twice the OMI sampling frequency. Figures 5 and 6 show the OMI 
ITFs and decompositions for the NO 2 fitting wavelength range (FWHM = 0.639 nm, 3.0 samples 
per FWHM). Note the asymmetry in the measured OMI ITFs and their decompositions. 

It is clear from all three examples that higher sampling by a factor of 2 greatly decreases 
the undersampled portion of the ITF, almost eliminating undersampling. Since for GOME this 
corresponds to 2.8 samples per FWHM, this demonstrates that information aliased from the 
spatial frequency band oo, nax < co < 2a) max represents the most important source of undersampling 
in this spectral region. 

5. Binning into Detector Pixels 

GOME and SCIAMACHY use Reticon-S RL-1024 SRU-type linear diode array detectors, which 
are designed to have the sensitivity profiles shown in Figure 7 for the spectral resolution and 
sampling of GOME. These detectors have 1024 photodiode elements spaced at 25 pm. The 
center 13 pm of each element is n-doped, with the bulk of the material p-doped, to generate the 
response profile shown (the newer Reticon-L detectors are also spaced at 25 pm, but with 19 pm 
n-doped regions and 6 pm interdiode gaps). 

Characterization of the actual pixel responses must take into account the fact that the 
input light source (mostly Fraunhofer spectrum, convolved with the ILS) varies significantly 
across the detector pixel profile. The final response is the integral across the pixel profile of the 
convolution of the input spectrum with the ILS: 


9 


( 13 ) 


**a)= 1 s(A*)®nA-w, 

prttftlr 


where R pix is the pixel response, S is the input spectrum, IP is the ILS, and ® denotes convolution. 
In order to calculate the response, it is first necessary to obtain the ILS by deconvolution, since 
the measured 1TF is the convolution of the ILS with the pixel response profile. This is 
accomplished using the Jansson method, as described in Ref. 14, for GOME, and a nonlinear 
least squares fitting to a parameterized profile shape for OMI. The result for the GOME slit 
function is shown in Figure 8. The deconvolved slit is the best match to the ILS that could be 
obtained in the iterative deconvolution process. The difference between the initial and 
reconvolved slit is an indicator of the goodness-of-fit for this procedure. The response functions 
for the OMI CCD detectors are Gaussian with a full-width at half maximum of 25 pm and 
separated by 22.5 pm (M. Dobber and R. Dirksen, private communication, 2003). Figures 9 and 
10 show these along with the decomposition to determine the ILS for the BrO and NO 2 fitting 
regions of OMI. 

6. Undersampling Correction 

6. 1 Wavelength calibration issues 

Satellite radiances and irradiances can be calibrated in wavelength with high absolute accuracy 
(typically <0.0004 nm for the GOME BrO fitting region) using cross-correlation to the solar 
spectrum of Ref. 3, with simultaneous fitting of the ITF. 212 Ref. 2 suggested that the apparent 
wavelength shift between GOME radiances and irradiances was due to instrumental effects, 
influencing the way the detectors are illuminated in the different measuring geometries. Ref. 4 
stated that it is due to the Doppler effect, since irradiances are obtained during the portion of the 
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orbit when the satellite is moving toward the sun. For ERS-2 orbit 81003031, the maximum 
velocity toward the sun is 7.46 km s ', while the average wavelength shift (irradiance - radiance) 
for the BrO fitting window is equivalent to 7.865±0.025 km s' 1 (0.0092310.00003 nm). The 
satellite Doppler shift contributes 6.89 km s' 1 (0.0081 nm), rather than the full 7.46 km s' 1 , given 
the GOME solar measurement procedure. 4 There is an additional component of 0.50 km s' 1 at 
this season from the ellipticity of the Earth’s orbit, for a total Doppler shift of 7.39 km s' 1 , 
leaving a significant instrument and/or spectral component (see the discussion in Ref. 15 for 
possible spectral “tilt” contributions). For the other spectrometers considered here, the relative 
contributions are not yet determined. For OMI in particular, since its undersampling is calculated 
here, a relative shift equal to the full Doppler shift is assumed; this will be modified when flight 
data become available. 

6.2 Undersampling calculations 

Calculation of the undersampling correction for each case is accomplished by convolving the 
high spectral resolution (0.01 nm) solar reference spectrum of Ref. 3 with the ITF determined for 
the instrument and wavelength region, and differencing fully-sampled and undersampled 
representations of this convolved solar spectrum at the sampling grid of the satellite radiances: 

1. Convolve the high resolution solar reference spectrum Eref with the satellite instrument 
ITF to create a lower spectral resolution, but highly oversampled, solar reference 
spectrum, E over : 

E ow M)= J E ref (A') ® ITF (A - A')dA . (14) 

profile 

The portion of E re f used in fitting BrO for GOME is shown in Figure 11, top panel; E over 
for this spectral region is shown in Figure 11, middle panel. 
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2. Determine from direct cross-correlation (for GOME) or estimate (for OMI) the 
wavelength grids for irradiance and radiance spectra g irr and g ra d- 

3. Sample E mer at the wavelength grid g irr , to give E Irr , and at the grid g rarf , to give E raii , 
using cubic spline interpolation 16 to determine values at each exact grid point. These are 
now undersampled representations of the solar reference spectrum, although each is 
correct at the points on its sampling grid. Eirr is shown in Figure 1 1, bottom panel; E ra d is 
virtually indistinguishable when over-plotted. 

4. Resample E irr to the wavelength grid g ra d , using cubic spline interpolation, giving E rad . 

5. The undersampling correction C„, in optical thickness units, is the difference between 
Erad and E! ra d, normalized to the average of E ra d over the fitting window: 

c.w-s-<a L-^ w . us) 

Erad 

C u corresponds to differencing fully-sampled and undersampled representations of the convolved 
solar spectrum since only E irr is resampled, thus inducing undersampling error on the radiance 
wavelength grid g rad . This undersampling correction assumes that, to first order, the radiance 
spectrum consists mostly of back scattered Fraunhofer structure. Higher-order corrections could 
be made to account for atmospheric absorption and the Ring effect. Figure 12 shows the 
residuals from fitting for BrO (without the use of an undersampling spectrum as a basis function) 
in ERS-2 orbit 81003031 for a measurement pixel midway through the orbit (pixel 800) and as 
an orbit average (top panel); the bottom panel shows the undersampling C u calculated here for 
the average relative wavelength shift between irradiance and radiance fitted for this orbit (0.0092 
nm). The calculated undersampling correction accounts for more than 90% of the fitting 
residuals. Figure 13, top panel, shows the undersampling calculated for the OMI baseline BrO 
fitting window, assuming a relative radiance-irradiance shift corresponding to the full Doppler 
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shift for the Aura orbit (7.5 km s' 1 ). The bottom panel shows the absorption for a typical BrO 
slant column density of lxlO 14 cm' 2 . The units (optical thickness and transmission) are 
equivalent in scale for these small interferences and absorptions. The undersampling correction 
is roughly an order of magnitude smaller than for GOME, but is still larger than the BrO 
absorption and needs to be carefully included in the fitting of the satellite data. Figure 14, top 
panel, shows the undersampling calculated for the OMI baseline NO 2 fitting window, assuming a 
relative radiance-irradiance shift corresponding to the full Doppler shift for the Aura orbit. The 
bottom panel shows the absorption for a typical NO 2 North American summer slant column 
density of 1.25xl0 16 cm' 2 . The contribution is less severe but, at ca. 20% of the nominal NO 2 
signal, it is still by no means negligible. Undersampling correction will very likely need to be 
included in order to derive meaningful tropospheric NO 2 abundances from OMI measurements. 
Under heavily polluted conditions over North America in summertime, for example, the 
tropospheric contribution to the NO 2 slant column usually does not exceed 40%. 6 Effective 
monitoring of moderate pollution requires correction to substantially better than that level. 

An attempt was made to improve the undersampling correction for GOME by calculation of 
the undersampling spectrum for the deconvolved ILS, followed by convolution of the resulting 
spectrum with the Reticon response function. The result of this procedure was expected to be an 
improved undersampling correction which would correspond even more closely to the GOME 
fitting residual, as shown in Figure 12, top panel. The actual result was substantially worse. The 
magnitude was approximately correct, but the spectral details were not. We think that the reason 
is that the actual detector response function deviates considerably from the trapezoidal shape 
shown in Figure 7. For the present, at least, the best undersampling corrections continue to be 
those determined using the entire ITF. These are shown in Figures 12, 13, and 14. When OMI 
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flight data become available, this procedure will be attempted to see whether it provides 
improved undersampling correction for this case. 

7. Discussion and Conclusions 

We provide a method that may be used to simply and reliably estimate the degree of 
undersampling an instrument configuration will have. For instruments where the ITF may be 
approximated as Gaussian (GOME), we show that sampling at 2.8 pixels/FWHM will almost 
completely eliminate undersampling. For more complex ITFs (OMI), significant undersampling 
can exist at 3.0 pixels/FWHM. 

Undersampling for ground-based, zenith-sky spectrometers has been discussed in Ref. 17 
in the context of measuring O 3 , NO 2 , and NO 3 . On the basis of numerical experiments, it is 
recommend to use sampling ratios between 4.5 and 6.5 pixels/FWHM, for Gaussian ITFs, to 
avoid undersampling. This is consistent with our finding that for OMI NO 2 measurements 
undersampling will be significant, and will highly impact tropospheric NO 2 measurements, at 3.0 
pixels/FWHM, but that at 6.0 pixels/FWHM it becomes negligible. The improvement in 
undersampling of OMI over that of GOME is due to the higher sampling rates (2.8 and 3.0 
samples per FWHM for the OMI UV and visible examples versus 1.4 for the GOME UV 
example). Asymmetric ITFs such as those of OMI and perhaps other imaging spectrometers may 
require higher sampling ratios than symmetric ITFs. 

The previously developed undersampling correction 2 4 is now commonly used in GOME 
scientific analyses and has been implemented operationally for SCIAMACHY. It has been 
demonstrated here that correction will be required for OMI BrO and tropospheric NO 2 
measurements, at least. Complete implementation will not be possible until the ITFs and 
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irradiance-radiance wavelength shifts are characterized in flight versus CCD row. A similar 
conclusion almost certainly applies to the OMPS instruments. 

Further refinement of the undersampling correction for OMI, to include full averaging of 
the solar spectrum over the ILS and convolution with the pixel response is currently in progress, 
to be ready for application to flight spectra. For GOME, this further refinement has been shown 
not to be an improvement, likely due to incorrect characterization of the Reticon response 
function. Convolution with the ITF followed by the sampling procedure described above in 
Undersampling calculations provides very accurate correction. 

Higher order corrections for instrumental effects, such as uneven sampling in wavelength 
space, and spectroscopic effects, such as atmospheric absorption and the Ring effect, may 
eventually be developed if analysis of flight data indicates that they are warranted. 
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Figure Captions 

Fig. 1 GOME Gaussian instrument transfer function in the spectral region used for bromine 
monoxide determinations, and its decomposition into Nyquist-sampled and undersampled 
components. 

Fig. 2 GOME Gaussian instrument transfer function and the Nyquist-sampled and undersampled 
portions for die hypothetical case where the slit function is sampled to twice the GOME 
spatial frequency. 

Fig. 3 OM! instrument transfer function for the BrO fitting region and the Nyquist-sampled and 
undersampled portions. 

Fig. 4 OMI instrument transfer function for the BrO fitting region and the Nyquist-sampled and 
undersampled portions for the hypothetical case where the slit function is sampled to twice 
the OMI spatial frequency. 

Fig. 5 OMI instrument transfer function for the NO 2 fitting region and the Nyquist-sampled and 
undersampled portions. 

Fig. 6 OMI instrument transfer function for the NO 2 fitting region and the Nyquist-sampled and 
undersampled portions for the hypothetical case where the slit function is sampled to twice 
the OMI spatial frequency. 

Fig. 7 Response of three Reticon-S detector pixels as a function of location in GOME detector 
channel 2. 

Fig. 8 Deconvolution of the GOME 1TF from the pixel response function to determine the ILS in 
the BrO fitting region. 

Fig. 9 Deconvolution of the OMI ITF from the pixel response function to determine the ILS in 
the BrO fitting region. 
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Fig. 10 Deconvolution of the OMI ITF from the pixel response function to determine the 1LS in 
the NO 2 fitting region. 

Fig. 1 1 The portion of the high resolution solar reference spectrum E re f used in fitting BrO for 
GOME (top); E re f convolved with the GOME ITF to create the lower spectral resolution, but 
highly oversampled, solar reference spectrum, E over (middle); E over sampled at the wavelength 
grid of the GOME irradiance (E irr ) (bottom). E irr and E ra d (not plotted here; see text) are 
undersampled representations of the solar reference spectrum. 

Fig. 12 Residuals from fitting GOME spectra for BrO in ERS-2 orbit 81003031, for a single 
spectrum and as an average for all spectra in the orbit (top); Synthetic undersampling, C u , 
calculated here for the relative wavelength shift between the GOME irradiance and the 
GOME radiance (bottom). 

Fig. 13 Synthetic OMI undersampling, C u , in the BrO fitting region, calculated here for a relative 
wavelength shift of the irradiance and radiance corresponding to 7.5 km s' 1 (top); absorption 
for a typical BrO slant column density of lxlO 14 cm' 2 (bottom). 

Fig. 14 Synthetic OMI undersampling, C u , in the NO 2 fitting region, calculated here for a relative 
wavelength shift of the irradiance and radiance corresponding to 7.5 km s' 1 (top); absorption 
for a typical N0 2 slant column density of 1.25x10 16 cm' 2 (bottom). 
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